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Measurements of the production of forward high energy 7r° mesons from transversely polarized 
proton collisions at y^i = 200 GeV are reported. The cross section is generally consistent with 
next-to- leading order perturbative QCD calculations. The analyzing power is small at if below 
about 0.3, and becomes positive and large at higher xf, similar to the trend in data at \fs < 20 
GeV. The analyzing power is in qualitative agreement with perturbative QCD model expectations. 
This is the first significant spin result seen for particles produced with pr > 1 GeV/c at a polarized 
proton collider. 

PACS numbers: 13.85.Ni, 13.88+e, 12.38.Qk 



An early qualitative expectation from perturbative 
Quantum Chromodynamics (pQCD) was that the chi- 
ral properties of the theory would make the analyzing 
power for inclusive particle production be very small 



The analyzing power (An) is the azimuthal asym- 
metry in particle yields from a transversely polarized 
beam incident on an unpolarized target. Earlier exper- 
iments studied polarized proton collisions (p-\ + p) at 
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center-of-mass energies yfs < 20 GeV and measured An 
for pion production at moderate transverse momentum 
(0.5 < pt < 2.0 GeV/c). Contrary to the naive expecta- 
tion, Am was found to be 20 — 40% for pions produced at 
large values of Feynman-x [xp = 2p£/Vs,wliere pl is 
the longitudinal momentum of the pion) |2|,|3j, |4], |5| . Sim- 
ilarly, elastic proton 6] and recent semi-inclusive deep- 
inelastic lepton scattering experiments 0, 0] have re- 
ported transverse single-spin asymmetries which differ 
significantly from zero. These results have sparked sub- 
stantial theoretical activity to understand transverse spin 
effects within the framework of pQCD 0. 

Pcrturbative QCD calculations of pion production in- 
volve the convolution of parton distribution and fragmen- 
tation functions with a hard partonic interaction. The 
reliability of calculations in the pQCD framework is ex- 
pected to increase with px- In this framework, forward tt 
production in p + p collisions is dominated by scattering 
of a valence quark in one proton from a soft gluon in the 
other. At large pseudorapiditics (77) and y/s < 20 GeV, 
there may be significant contributions to particle pro- 
duction from soft hadronic processes collectively known 
as beam fragmentation. At a collider, y/s is significantly 
larger, leading to the expectation that the origin of for- 
ward pions will shift towards collisions of the partonic 
constituents of the proton, consistent with the PYTHIA 
Monte Carlo generator |l0j . Measurements of the cross 
section for forward pion production are important to es- 
tablish that pQCD is a suitable framework for treating 
polarization observables in these kinematics. 

Different mechanisms have been identified in the 
pQCD framework by which one might exp ect transverse 
spin effects El El El El El OS E3, all of which 
may contribute to some degree. With only data at 
yfs < 20 GeV for comparison, these models are not well 
constrained. Despite this, the models have been extrapo- 
lated by an order of magnitude in -y/s and approximately 
a factor of 2 in pt, and all predict that sizable transverse 
spin effects will persist at y/s = 200 GeV. This Letter 
addresses the question if An is sizable at y/s — 200 GeV, 
as predicted by these models. We present measurements 
of the cross section and An for the production of forward 
7T° mesons having px > I GeV/c from p^ + p collisions 
at y/s = 200 GeV. 

Data were collected by the STAR experiment (Solenoid 
Tracker at RHIC) at the Brookhaven National Labora- 
tory Relativistic Heavy Ion Collider (RHIC) in January 
2002. RHIC is the first polarized proton collider. Po- 
larization is produced by optical pumping of an atomic- 
beam source E3 and is partially preserved through an 
accelerator complex to reach RHIC E3- ^ n RHIC, a pair 
of helical dipole magnets in each ring serves as the first 
use of full "Siberian snakes" [3] in a high-energy accel- 
erator to preserve polarization during beam acceleration. 
The stable spin axis of the RHIC rings is vertical. Beam 
bunches crossed the STAR interaction region (IR) every 



213 nsec. The polarization direction alternated between 
up and down for successive bunches of one beam and after 
every two bunches of the other beam. Data were sorted 
according to the spin direction of the beam correspond- 
ing to positive xf pion production. Summing all bunches 
in the other beam resulted in negligible remnant polar- 
ization. Typical luminosities were 10 30 cm _2 sec _1 , and 
the integrated luminosity was 150 nb _1 for these data. 

The average beam polarization for each fill, Pbeam, was 
determined using a Coulomb-Nuclear Interference (CNI) 
polarimeter located in RHIC HIH^]. At 24.3 GeV, the 
RHIC injection energy, the analyzing pow er of the CNI 
reaction is A% NI = 0.0133 ± 0.0041 [H 0, and can 
be used to deduce the absolute polarization of the pro- 
ton beam. However, at 100 GeV, the beam energy used 
for RHIC collisions, A N NI has not yet been measured. 
The CNI asymmetries measured at injection and colli- 
sion energies were nearly equal for many fills. Since the 
beam acceleration process is unlikely to increase Pbeam, 
this suggests that A N NI at 100 GeV is no smaller than 
at 24.3 GeV. For the present analysis, we assume there 
is no change in A N NI between these two energies, giving 
an average value of (Pbeam) — 0.16. 

A prototype forward 7r° detector (PFPD) was installed 
at STAR 750 cm from the IR to identify ir° mesons. At 
this time, STAR does not have the capability to iden- 
tify large rapidity charged pions. The PFPD consisted 
of a Pb-scintillator sampling calorimeter [25j , placed with 
its edge ~ 30 cm left of the oncoming polarized proton 
beam (beam-left). The PFPD was 21 radiation lengths 
deep and subdivided into 4x3 towers. To measure the 
transverse profiles of photon showers, the PFPD had a 
shower-maximum detector (SMD) approximately 5 radi- 
ation lengths deep, comprising two orthogonal layers of 
100 x 60 scintillator strips spaced at 0.37 cm. To ad- 
dress systematic errors associated with measuring left- 
right asymmetries with a single arm detector, an array of 
Pb-glass detectors with no SMD was placed to the right of 
the oncoming polarized proton beam (beam-right). Sim- 
ilar arrays were placed above and below the vertically 
polarized beam, where no asymmetries are expected. 

The luminosity was measured at STAR using beam- 
beam counters (BBC) [2^ composed of segmented scin- 
tillator annuli mounted around the beam at longitudi- 
nal positions z — ±370 cm, spanning 3.3 < < 5.0. 
Proton collision events were identified by requiring the 
coincidence of at least one BBC segment fore and aft 
of the IR. Absolute luminosity was determined by mea- 
suring the transverse size of the colliding beams and 
the number of protons colliding at STAR. The cross 
section measured for the BBC coincidence condition is 
26.1 ± 0.2(stat.) ± 1.8(svst.) mb [2*3, consistent with sim- 
ulation [ljj E3|. The BBC observes 87 ± 8% of the in- 
elastic, non-singly diffractive cross section. 

All forward calorimeters were read out when the energy 
deposited in any one calorimeter was greater than that 
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from an electron of ~ 15 GeV. The BBC coincidence 
requirement was imposed to select p + p collisions. 
The asymmetry measured at beam-left is 



P 



beam-^N 



N 4 



(1) 



The number of tt° mesons detected when the beam spin 
vector was oriented up (down) is N+r_y The spin- 
dependent relative luminosity (R = C+/C- ~ 1.15) was 
measured with the BBC. Background contributions to R 
were reduced by increasing the coincidence requirements 
to at least two BBC segments on each side of STAR. 
The systematic errors on R, primarily coming from the 
change in R when the background is corrected, are of the 
order of 10~ 3 and are a factor of 10 to 20 smaller 
than PbeamAN measured with the PFPD. 

Neutral pions are reconstructed utilizing the formula 

M 77 = En^Jl - sin(</> 77 /2) w E tot J\ — z* d~ n /2z vtx 

using events with at least two clusters in the SMD. The 
energy of the leading 7r°, E w , is taken to be the total 
energy deposited in all of the towers, E tot . The opening 
angle between the two photons, 77 , is determined by 
z v txi the distance between the collision vertex and the 
PFPD, and the separation of the photons at the detector, 
<i 77 . Both rf 77 and z 1 — \E 1 \ — E 1 2\/(E 1 i + E-/2) are 
measured by an analysis of the energy deposited in the 
strips of the SMD. The value of d 77 is determined from 
the fitted centroids of the peaks, while z~ i is derived from 
the ratio of the fitted areas under the peaks. A fiducial 
volume is defined by requiring the SMD peaks to be more 
than 12 strips from the detector edge. Figure^shows the 
M 77 spectra for two energy bins. The mass resolution is 
20 MeV/c 2 (RMS) for 15 < E tot < 80 GeV, limited by 
the measurement of 77 . The centroid of the ir° peak is 
used to determine the calibration for each tower for each 
fill to an accuracy of the order of 1%. The calibration is 
found to have negligible dependence on energy or spin- 
state. 

The 7T° detection efficiency is determined in a matrix 
of E n and 77 from a Monte Carlo (MC) simulation of p+p 
collisions and the detector response [2^. The open 
histograms in Fig. ^ are MC events which undergo the 
same reconstruction and selection as the data. The MC 
matches the data well for several variables, including px, 
E to t, and 77. The tt° detection efficiency is dominated by 
the geometrical acceptance of the calorimeter. 

The 7T° sample is distorted by coincident particles 
from the jets containing them. The PFPD is about one 
hadronic interaction length deep. When two photons 
from 7T° decay overlap with other particles, the PFPD 
response to the other particles tends to increase E to t rel- 
ative to En and broaden the 77 resolution. This results 
in a broad M 77 distribution peaked at a value larger 
than M„. The average value of E tot is about 3 GeV 
larger than E n , independent of E w . MC events with 
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FIG. 1: Uncorrected spectra of the diphoton invariant mass 
in two energy bins. The points are data with statistical 
errors. The open histograms are reconstructed simulation 
events, normalized to equal area. The hatched histograms 
are used to correct the cross section. 



\E tot —E % \ > 2 GeV are shown as the hatched histograms 
in Fig.^ Events with only one photon from 7r° decay plus 
other particles exist predominantly at small M 77 , and 
are suppressed by requiring z 7 < 0.3. The ^-dependent 
systematic error in the cross section is about 20%, domi- 
nated by the jet correction. The MC simulation includes 
7T° events from forward jets. The uncertainty includes 
the difference when these effects are explicitly corrected 
in both the data and the simulation, and in neither. 

Non-collision background is suppressed to the level of 
1% by requiring the coincidence from the BBC in the 
offline analysis. Following our simulations, the cross sec- 
tion is corrected by 10% to account for the bias intro- 
duced by the BBC coincidence condition. Hadronic back- 
ground comprising events with no leading 7r° in the accep- 
tance of the calorimeter is predominantly at small M 77 , 
and is reduced by constraining z 1 . The hadronic back- 
ground amounts to about 2% of the yield underneath the 
7T° peak at 0.09 < M 77 < 0.22 GeV/c 2 . 

The inclusive ir° production cross section for 30 < 
En < 55 GeV in 5 GeV bins is presented in Fig. |21 Data 
with 3.4 < r\ < 4.0 were selected, giving (77) = 3.8 in- 
dependent of E^; in this range the detector efficiency 
is well understood. The dominant contributions to the 
normalization error come from knowledge of the absolute 
transverse position of the detector (10%), the absolute lu- 
minosity determination (8%), and the model dependence 
of the BBC efficiency (8%). The data are plotted at the 
average E„ of the bin. 

The curves on the plo t are next-to-leading order (NLO) 
pQCD calculations [29| evaluated at 77 = 3.8, using the 
CTEQ6M |32j parton distribution functions and equal 



5 



> 
(5 



10 



£2 



Q. 

73 

e 

LU 



10 





Data 3.4<r|<4.0 




• 7i° mesons «r|>=3.8) 


- v\ 


NLO pQCD calc. 




- KKP F.F. (ri=3.8) 




- Kretzer F.F. (ri=3.8) 


r 

- Normalization 




Uncertainty = 


17%\ \ 


:<p T >= 1.5 1.7 1.9 2.1 '••?.2\GeV/c 



25 30 35 40 45 50 55 60 65 

(GeV) 

FIG. 2: Inclusive 7r° production cross section versus leading 
7T° energy (E n ). The average transverse momentum ((pt)) 
is correlated with E n , as the PFPD was at a fixed pseudo- 
rapidity (j]). The inner error bars are statistical, and are 
smaller than the symbols for most points. The outer error 
bars combine these with the iJ^-dependent systematic errors. 
The curves are NLO pQCD calculations evaluated at r\ = 3.8 




FIG. 3: Analyzing powers versus Feynman-a; (xf). The aver- 
age transverse momentum {{pt}) is correlated with xf- The 
solid points are for identified tt° mesons. The open points are 
for the total energy (Etot), shifted by xf + 0.01. The inner er- 
ror bars are statistical, and the outer combine these with the 
point-to-point systematic errors. The curves are from pQCD 
models evaluated at pr = 1.5 GeV/c O EH E3 The An 
values are proportional to A N NI , assumed to be 0.013 at 100 
GeV. 



renormalization and factorization scales of pt- The NLO 
pQCD calculations are in general consistent with the 
data, in contrast to midrapidity ir° data at lower y/s 
The solid line uses the "Kniehl-Kramer-Potter" 
(KKP) set of fragmentation functions (F.F.) [3fJ, while 
the dashed line uses the "Kretzer" set [2l|. The differ- 
ence between the two reflects uncertainties in the F. F. 
at these kinematics. At the chosen scale, KKP tends to 
agree with the data better than Kretzer, consistent with 
midrapidity tt° data at s/s = 200 GeV 

The analyzing power is presented in Fig. EI plotted ver- 
sus 2 (Etot)/y/s ~ xp. The solid points are for tt° mesons 
from 3.3 < r) < 4.1 and 0.07 < M 77 < 0.3 GeV/c 2 , 
with a^-dependent constraints on z 7 to minimize back- 
ground. The open points are based solely on E to t in the 
PFPD without SMD analysis: neither fiducial volume 
constraints nor 7r° identification. The agreement between 
the solid and open points indicates An is not sensitive to 
the analysis used to identify 7r° mesons. This is consis- 
tent with simulations showing that 95% of events with at 
least 25 GeV deposited in the PFPD come from photons, 
95% of which are daughters from tt° decay. The An seen 
at beam-right with the Pb-glass array is similar to that 
seen at beam-left with the PFPD, while An for the Pb- 
glass above and below the beam is consistent with zero, 
as expected. The largest x_F-dependent systematic error 
arises from comparison of the beam-left and beam-right 
data. The average An(xf) is computed using both, and 
an uncertainty is assigned to bring An(xf) (shown in 



Fig. |2J) within 1 standard deviation of the average. 

The curves on the plot are predictions from the pQCD 
models, fitted to data at y/s — 20 GeV, extrapolated 
to y/s = 200 GeV and evaluated at pr — 1.5 GeV/c 
[II El El 03 One model attributes single-spin effects 
to the convolution of the transversity distribution func- 
tion with a spin-dependent Collins fragmentation func- 
tion [l^j . The Sivers model adds explicit spin-dependent 
Ut dependence to the parton distribution functions |l5l |. 
Other models ascribe the effects to twist-3 parton corre- 
lations in the initial or final state 0, ^3 • The data are 
qualitatively consistent with all of these predictions. 

The trend of An at lower y/s is to increase from zero 
beginning at a value of xf which depends on y/s |5j . The 
significance of the increase for these data is 4.7 o (includ- 
ing statistical and point-to-point systematic errors) from 
a linear fit to the open circles in Fig.|21for xf > 0.27, with 
X 2 = 0.9 for 3 degrees of freedom. This is the first sig- 
nificant spin result seen for particles with pt > 1 GeV/c 
at a polarized proton collider. 

In summary, high energy 7r° mesons have been ob- 
served at forward angles from p-\ + p collisions at i/s = 
200 GeV. The differential cross section is, in general, 
consistent with NLO pQCD calculations. The analyz- 
ing power is small at xf below about 0.3, and becomes 
positive and large at higher xf, similar to the trend ob- 
served in fixed-target data at y/s < 20 GeV. The analyz- 
ing power at y/s — 200 GeV is in qualitative agreement 
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with pQCD model predictions. Higher precision mea- 
surements of An as a function of both xp and px may 
help to differentiate among the models. Future measure- 
ments may attempt to determine the Collins fragmenta- 
tion function in pj +p collisions, as well as to look at jet 
production and Drell-Yan scattering to isolate potential 
contributions to transverse spin effects. 
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